AD-A037  845 


UNCLASSIFIED 

of  | 

40A037845 


CINCINNATI  UN IV  OHIO  DEPT  OF  ENGINEERING  ANALYSIS  F/G  6/5 

UCIN  VEHICLE-OCCUPANT/CRASH-VICTIM  SIMULATION  MODELf(U) 

JUN  76  R L HUSTONf  C E PASSERELLOf  M V HARLOW  N00014-72-A-0627-0002 


OOC  file  COPY  __  flOA 0 37 8 45 


UCIN  VEHICLE-OCCUPANT/CRASH-VICTIM  SIMULATION  MODEL v 
L-..  J 


Report  Date;  I June  197_6 

R.  L. /Huston,  £•  E. /^asserello,' *»d  M.  W.  barlow 

University  of  Cincinnati 
Cincinnati,  Ohio  45221 

Til ( 

Contract  No. I N0001if-72-A-OO27-OOO2 


-L  /"■  '/<i\ 


/ i 


% \ 


INTRODUCTION 


Recently  there  have  appeared  a number  of  finite-segment  computer 
models  of  the  human  body  designed  primarily  to  conduct  crash-victim 
analyses.  ; Indeed,  there  are  currently  as  many  as  10  distinct  vehicle- 
occupant/crash-victim  models  available  [1],  Perhaps  the  first  of 
these  was  a two-dimensional  computer  code  developed  originally  by 
McHenry  [2]  in  1963  and  later  refined  in  1966  and  1968.  More  recent 
two-dimensional  models  include  those  developed  by  Segal  [3],  Danforth 
and  Randall  [4],  Robbins,  Bowman  and  Bennett  [5],  Glancy  and  Larsen 
[6],  and  Karnes,  Tocher  and  Twigg  [7].  The  basic  difference  in  these 
codes  is  the  variety  of  input-output  options  available.  The  three- 
dimensional  models  include  those  developed  by  Robbins  [8],  Robbins, 
Bennett  and  Bowman  [9],  Young  [10],  Furusho  and  Yokoya  [11],  Bartz 
112],  Fleck,  Butler  and  Vogel  [13,  14],  and  Huston,  Passerello, 

Harlow,  Hessel  and  Winget  [15-18].  There  is  more  diversity  among 
these  codes  because  of  their  greater  complexity.  A summary  of  the 
various  characteristics  of  these  codes,  their  capabilities,  and  their 
formulation  may  be  found  in  [1]. 

The  objective  of  this  paper  is  to  present  a description  of  the 

UCIN  code  with  particular  emphasis  upon  its  options,  capabilities, 
formulation  and  input-output  procedures.  The  balance  of  the  paper 
is  divided  into  five  parts  with  the  first  of  those  providing  a 
description  of  the  model  and  its  capabilities.  This  is  followed  in 
the  next  three  parts  by  listing  of  the  input-output  procedures  and  the 
final  part  provides  a summary  of  the  vital  features  of  the  code. 


MODEL  DESCRIPTION  AND  CAPABILITIES 

The  model  consists  of  12  rigid  bodies  representing  the  human 
limbs  together  with  a vehicle  cockpit  as  shown  in  Figure  1.  The 
twelve  bodies  of  the  model  are  connected  together  with  ball-and- 
socket  joints  as  shown  in  Figure  2. 
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Figure  2.  Side  View  of  the  Model  in  Sitting  Configuration 

Forty-eight  variables  are  required  to  describe  the  position  and 
orientation  of  the  model.  These  are: 


4’X5’X6 

7,X8,X9 


X10,X11,X12 


X13,X14,X15 


X16*X17,X18 


X19*X20,X21 


X22,X23’X24 


position  of  the  vehicle  relative 
to  an  Inertial  frame 

orientation  of  the  vehicle 
relative  to  an  inertial  frame 

position  of  a reference  point 
in  B2,  the  lower  torso  relative  to 
the  origin  of  the  vehicle  frame 

orientation  of  B2,  the  lower  torso, 
relative  to  the  vehicle  frame 

orientation  of  B3,  the  middel  torso, 
relative  to  B2,  the  lower  torso 

orientation  of  B4,  the  upper  torso, 
relative  to  B3,  the  middel  torso 

orientation  of  B^,  the  upper  left  arm, 
relative  to  B5,  the  upper  left  arm 

orientation  of  Bg,  the  lower  left  arm 
relative  to  B5,  the  upper  left  arm 


X25,X26,X27 

orientation  of  By,  the  upper  right 
arm,  relative  to  B4,  the  upper  torso 

X28,X29,X30 

orientation  of  Bg,  the  lower  right 
arm,  relative  to  B7,  the  upper  right  arm 

X31’X32,X33 

orientation  of  Bg,  the  head,  relative  to 
B4,  the  upper  torso 

X34’X35,X36 

orientation  of  Bj«,  the  upper  left  leg. 

relative  to  B2,  tne  lower  torso 

X37,X38,X39 

orientation  of  B^,  the  lower  left  leg, 
relative  to  B^g,  the  upper  left  leg 

x40*x41‘x42 

orientation  of  B-j^,  the  upper  right  leg, 
relative  to  B2,  the  lower  torso 

x43*x44»x45 

orientation  of  B^_,  the  lower  right  leg, 
relative  to  B^2»  the  upper  right  leg 

X46,X47,X48 

position  of  a reference  point  in  B9,  the 
head,  relative  to  an  attach  point  in  B4, 
the  upper  torso 

All  of  these  variables,  except  for  the  position  variables,  x^, 

X2,  X3,  X7,  xg,  X9,  X43,  X47  and  x^g,  are  orientations  angles  generated 
by  successive  rotation  of  adjacent  bodies  relative  to  each  other.  The 
model  thus  uses  "relative  coordinates"  (ie.  orientation  between  bodies) 
as  opposed  to  "absolute  coordinates"  (orientation  in  space)  to  define 
its  configuration.  These  orientation  angles  are  in  turn  defined  as 
follows:  Imagine  two  adjacent  bodies  of  the  system  oriented  relative 

to  each  other  so  that  coordinate  axes  (X,Y,Z)  unbedded  in  each  are 
respectively  aligned.  Next,  imagine  the  second  body  rotated  dextrally, 
relative  to  the  first  body,  through  an  angle  a,  about  its  X axis.  Then 
imagine  successive  dextral  rotations  of  the  second  body  about  its  Y and 
Z axes  through  angles  8 and  y.  This  brings  the  second  body  into  general 
orientation  relative  to  the  first  defined  by  the  relative  orientation 
angles  a,  8,  and  y.  The  orientation  variables,  listed  above  in  sets  of 
three,  are  respectively  a,  8,  and  y relative  orientation  angles. 

The  first  six  variables  listed  above  define  the  motion  of  the 
cockpit  or  vehicle  frame  relative  to  inertia  space.  These  variables 
must  be  specified  in  input  data.  Also,  variables  X22»  x24*  x28’  x30* 

X37,  X39,  X43,  and  X45  are  usually  specified  as  zero  to  simulate 
hinge  joints  at  the  elbows  and  knees.  The  remaining  31  variables  may 
be  either  specified  or  left  as  unknowns.  If  the  variables  are  specified 
(eg.  X36  * 0),  the  required  moment  needed  to  maintain  that  specification 
(eg,  Mjg)  is  determined. 

The  model  allows  for  the  arbitrary  specification  of  external  forces 
and  moments  on  each  of  its  bodies.  These  forces  and  moments  are 
represented  on  each  body  as  a single  force  passing  through  its  mass 
center,  together  with  a couple.  These  forces  and  the  moments  of  the 
couples  are  part  of  the  input  data. 


The  model's  initial  position  is  generally  in  an  erect  sitting 
position  as  shown  in  Figures  2 and  3.  In  this  poistion,  all  the  body 
coordinate  axes  and  the  vehicle  frame  are  aligned  and  thus  all  the 
orientation  angles  are  zero.  The  model  has  a seat  which  in  turn  is 
modelled  by  springs  as  shown  in  Figure  2.  Essentially,  there  are 
seven  springs  which  may  exert  forces  on  the  model  with  the  points 
of  contact  being  the  mass  centers  of  bodies  2,  3,  4,  9,  10  and  12. 
Viscous  damper  stops  are  used  to  provide  limited  seat  deflection. 

The  amount  of  seat  deflection  and  spring  constants  are  input  data. 


Figure  3.  The  Model  and  Seat  in  Initial  Configuration 

The  computer  model  provides  for  the  use  of  up  to  ten  restrain- 
ing belts  attached  at  arbitrary  points  between  the  cockpit  and  the 
bodies  of  the  model.  These  belts  are  modelled  as  springs  with  the 
spring  constants  and  attach  points  as  part  of  the  input  data. 

The  ball-and-socket  connection  joints  of  the  model  are  provided 
with  angle-stops,  modelled  by  one-way  dampers,  to  simulate  motion 
constraints  of  the  human  limbs.  The  angle  limits  and  the  damping 
coefficients  are  part  of  the  input  data. 

The  model  also  has  the  provision  of  a "stretching"  or  extending 
neck.  Neck  stretch  is  modelled  by  a spring  and  a damper  between 
attach  points  in  bodies  4 and  9.  The  spring  constant  and  the  damping 
coefficient  are  input  data. 

Finally,  the  computer  model  provides  for  the  use  of  12  intrusion 
surfaces  to  simulate  the  cockpit  or  vehicle  interior.  These  intrusion 
surfaces  are  planes  and  the  output  of  the  code  records  when  a body  or 
limb  of  the  model  collides  with  or  passes  through  an  intrusion  surface. 
The  Intrusion  surfaces  are  listed  as  follows: 


1) 

Left  windshield 

5) 

Lower 

left  i 

door 

9) 

Firewall 

2) 

Front  windshield 

6) 

Upper 

right 

door 

10) 

Top  dash 

3) 

Right  windshield 

7) 

Lower 

right 

door 

ID 

Front  dash 

4) 

Upper  left  door 

8) 

Roof 

12) 

Bottom  dash 

The  locations  and  inclination  of  these  intrusion  surfaces  are  also 
part  of  the  input  data. 


COMPUTER  CODE : INPUT/OUTPUT 


The  following  brief  paragraphs  provide  a general  description  of  the 
input  data  required  and  the  output  provided  by  the  computer  code.  The 
specific  format  used  is  described  in  detail  in  the  next  two  parts  of  the 
paper. 


The  input  consists  of  the  following: 

Physical  parameters.  These  are  the  masses,  inertia  dyadics,  mass  center 
positions  connection  point  positions,  and  orientation  angle  limits,  for 
each  of  the  12  bodies  of  the  human  model. 

Cockpit  geometry.  This  consists  of  a normal  vector  and  a location  point 
for  each  of  the  12  intrusion  surface  planes.  Also  the  floor  position  and 
supporting  spring  force  constant  are  part  of  the  cockpit  specifications. 

Cockpit  motion.  The  cockpit  displacement  and  rotation  relative  to  an 
inertia  frame  are  required  as  input  (x^,...,x^).  Typically,  it  is  desired 
to  express  this  in  the  form  of  linear  and  angular  acceleration  of  the 
cockpit.  The  computer  program  is  written  so  that  the  cockpit  accelera- 
tion components  may  be  "read  in"  by  simply  specifying  the  acceleration 
at  selected  time  intervals  of  the  acceleration  profile.  (This  is,  in 
effect,  a straight  line  approximation  to  an  acceleration  curve).  Six 
(three  translation  and  three  rotation)  such  acceleration  profiles  or 
curves  may  be  employed. 

Spring  and  damping  constants.  This  includes  the  seat,  restraining  belt, 
orientation  angle  constants,  and  neck  parameters.  Also,  the  attach 
points  of  the  restraining  belts  are  included. 

Initial  conditions.  This  includes  the  initial  values  of  the  unknown 
variables  and  their  derivatives.  Also  the  external  forces  and  moments 
(if  any)  which  are  applied  to  the  bodies  of  the  models  must  be  specified. 


The  output  consists  of  two  parts:  The  first  part  is  simply  an 

"echo"  or  copy  of  the  input  data.  The  second  part  contains  at  each 
output  time  the  following: 

1.  The  values  of  all  variables  and  their  first  derivatives. 

2.  The  joint  and  mass  center  positions  in  both  inertia  space 
and  relative  to  the  vehicle. 

3.  The  mass  center  velocities  and  accelerations. 

4.  The  moments  and  forces  associated  with  variables  which  are 
specified. 

5.  Restraining  belt  forces. 


SPECIFIC  INPUT  DATA 


Data  may  be  entered  into  the  program  either  by  following  the 
instructions  of  the  interrogating-terminal  version  of  the  program  or 
by  batch  (or  remote  job  entry).  The  following  paragraphs  provide  the 
details  of  the  required  input  card  images  for  batch  input.  This  is 
followed  by  a sample  listing  of  input  data. 

1)  Error  Set 

The  first  card  is  a blank  card.  (It  is  associated  with  a 
debugging  subroutine  of  the  program.  This  card  must  be 
included,  but  being  blank,  it  disengages  the  debugging 
routine) . 

2)  Mass  of  the  12  Bodies 

Next,  there  are  12  data  entries  on  3 cards  for  the  masses 
of  the  12  bodies.  They  are  read  sequentially  for  the  12 
bodies  in  a 5F10.9  format.  The  usual  units  are  slugs. 

3)  Inertia  Dvadics 

Following  this,  there  are  12  (3  x 3)  inertia  matrices  read 
sequentially,  one  for  each  body.  Each  matrix  requires 
three  cards.  Each  card  contains  a row  of  an  inertia  matrix. 
The  matrices  are  referenced  to  the  axis  system  as  shown  in 
Figure  A.  They  are  read  in  a 5F10.9  format.  The  usual 
units  are  slugs-in^.  These  are  a 36  cards  for  this  entry. 

A)  R Vectors 

There  are  12  of  these.  They  are  vectors  referenced  to  the 
body  co-ordinate  systems  which  locate  the  mass  center  of  a 
body  relative  to  the  connection  point  for  that  body  with  its 
adjacent  lower  numbered  body,  or  in  the  case  of  B2  relative 
to  the  chosen  references  point  in  B2.  (See  Figure  5).  They 
are  read  sequentially  one  vector  per  card  in  a 5F10.9  format. 
There  are  12  cards.  The  usual  units  are  inches. 

5)  Xl-Vectors 

There  are  eleven  of  these.  They  are  vectors  referenced  to 
body  co-ordinate  systems  which  locate  the  relative  positions 
of  connection  points  (See  Figure  5).  They  are  read 
sequentially  one  vector  per  card  in  a 5F10.9  format.  The 
usual  units  are  inches.  There  are  eleven  cards. 

6)  Cockpit  Geometry-Intrusion  Surfaces 

The  first  card  (the  first  card  in  this  sequence)  contains 
either  a zero  of  12  read  in  15  format.  If  a zero  is  read 
the  intrusion  surface  option  is  omitted,  and  no  additional 
cards  should  be  read  in  this  section.  If  a 12  is  read,  12 
additional  cards  should  be  read  (one  for  each  of  the  12 
intrusion  planes). 


Figure  5. 

v»efora  in  the  Model 


One  each  of  these  cards  six  values  are  read.  The  first 
three  values  are  the  coordinates  (X,Y,Z)  of  any  point 
in  the  intrusion  plane  relative  to  0,  the  cockpit  origin. 
The  last  three  values  are  the  components  of  the  outward 
directed  vector  normal  to  the  intrusion  plane.  The 
format  is  6F10.9.  The  12  intrusion  planes  are  to  be 
entered  in  the  following  order: 


Left  windshield 
Front  windshield 
Right  windshield 
Upper  left  door 
Lower  left  door 
Upper  right  door 


Lower  right  door 
Roof 

Firewall 
Top  dash 
Front  dash 
Bottom  dash 


(If  a number  less  than  12  but  greater  than  zero  is  read 
on  the  first  card,  then  the  intrusion  plane  geometry  for 
those  planes  not  read  will  default  to  the  YZ  plane). 


Magnetic  Tape  Card 

This  card  allows  the  user  to  write  all  the  output  data  on 
magnetic  tape  by  writing  the  tape  unit  number  in  15  format. 
If  a zero  is  read,  or  a blank  card  is  used,  no  magnetic 
tape  is  written. 

Heading  Card 

On  this  card  the  user  writes  anything  he  desires  as  a' 
heading . 


Floor  Spring  Constant  and  Distance 

This  is  a spring  constant  for  the  floor  and  the  vertical 
distance  from  the  reference  point  in  body  2 to  the  floor 
of  the  vehicle  when  the  model  is  in  the  reference  position. 
They  are  read  sequentially  on  a single  card  in  a 5F10.9 
format.  The  usual  units  are  lbs/in  and  inches. 


Acceleration  Input  for  the  Cockpit 

The  computer  program  has  been  written  so  that  the  cockpit 
acceleration  components  (three  translation  and  three 
rotation)  may  be  "read  in"  by  simply  specifying  the 
acceleration  at  selected  time  intervals  of  an  acceleration 
profile.  The  program  has  the  capability  of  including  up 
to  25  points  of  the  acceleration  curve. 


To  be  more  specific,  consider  a straight  line 
approximation  to  the  acceleration  curve.  The 
acceleration,  velocity  and  displacement  during 
the  i^  time  interval  are  then: 
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i+1  i 


ai+l”ai 

v * v + a (t-t  ) + (- ~)  -±- 
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d - di + V'-'i’ + ai— — + (i— — <r~ 


where  a^,  v^,  d^  and  t£  are  the  acceleration,  velocity, 
displacement  and  time  at  the  beginning  of  the  ic^  interval. 
The  kinematical  profile  of  the  cockpit  is  thus  determined 
when  the  a^  are  specified  and  when  v(and  dr,  the  initial 
velocity  and  displacement  are  given. 

The  computer  cards  are  coded  as  follows:  The  first  card 
contains  the  number  of  profiles  (up  to  6)  that  are  used 
in  1015  format. 

The  next  cards  contain  the  number  of  points  (up  to  25) 
and  the  variable  numbers  (i.e.  1 for  x,  2 for  y,  3 for  z, 

4 for  a,  5 for  6,  6 for  y)  of  the  respective  acceleration 
profile  curve  for  each  acceleration  profile  used.  They 
are  read  sequentially  in  a 1015  format. 

The  next  group  of  cards  contain  the  time  and  acceleration 
ordinate  for  each  point  of  the  respective  profile  curves 
which  are  used.  The  first  card  of  each  subgroup  (i.e. 
for  each  profile  used)  also  contains  the  initial  velocity 
and  displacement  for  that  profile.  The  data  on  these 
cards  are  in  5F10.9  format. 

11)  Neck  Spring  and  Damping  Constants 

These  constants  must  be  included  in  the  input  data  even 
if  the  neck  is  to  remain  unstretched  during  the  motion. 

(A  blank  card  may  be  used  if  the  neck  is  to  remain 
unstretched).  These  constants  ( spring  and  damping 
respectively)  are  entered  on  a single  card  in  a 5F10.9 
format.  The  units  are  usually  lb/in  and  lb. sec/in 
respectively. 


12)  Orientation  Angle  Limits 

There  are  66  of  these — a maximum  and  a minimum  angle 
displacement  for  each  of  the  33  internal  orientation 
angles.  For  example,  the  first  card  defines  the  limits 
on  the  orientation  of  B3  relative  to  B2  by  containing 
the  following  data: 


^°23  max’®23  max,Y23  max,a23  min’ ^23  max’^23  min^ 
or  (see  the  variable  list) 


^X13  max’X14  max’X15  max,x13  min,X14  min,X15  min^ 

Each  card  contains  6 data  entries  in  a 6G10.9  format  in 
the  above  order.  There  are  a total  of  11  cards.  The 
units  are  degrees. 

13)  Damping  Coefficient  for  the  Joints  (CDAKP1) 

This  is  the  value  of  the  damping  coefficient  to  be  used  for 
all  connecting  joints  of  the  human  body  model.  It  is 
entered  on  a single  card  in  a 5F10.9  format.  The  usual 
units  are  lb  sec/rad. 

14)  Damping  Coefficient  for  the  Joint  Limits  (CDAMP2) 

This  is  usually  a large  number  which  is  required  to  model  a 
joint  angle  displacement  unit.  It  is  read  on  a single  card 
in  a 5F10.9  format.  The  usual  units  are  lb  sec/rad. 

15)  Spring  Constants  for  the  Seat  (ESEAT[21) 

There  are  two  spring  constants:  one  for  the  seat  and  one 

for  the  back  of  the  seat.  ESEAT1  is  for  the  back  of  the 
seat  and  ESEAT2  is  for  the  seat.  They  are  read  sequentially 
on  a single  card  in  a 5F10.9  format.  The  usual  units  are 
lbs/in. 

16)  Displacement  Limits  for  the  Seat  (SEATM[2]) 


There  are  two  max  deflections:  one  for  the  back  of  the 
seat  and  one  for  the  bottom  of  the  seat.  They  are  read 
sequentially  on  a single  card,  in  a 5F10.9  format.  The 
usual  units  are  inches. 


switch 

5) 

Floor  constraints 

switch 

6) 

Impulses — (always  should  be 

switch 

7) 

Neck  forces — (usually  off) 

switch 

8) 

Intrusion  surfaces  (cockpit 

8 cards  are  read 

off  [zero]) 


geometry) 


20)  NOUT 


The  total  number  of  unknown  variables  (maximum  of  34)  is 
read  on  a single  card  in  15  format. 


21)  JFX1  Array 

JFX1  is  the  array  of  unknown  variables.  The  total  number 
of  these  is  NOUT.  The  unknown  variable  indices  (eg. 

. index  is  17)  are  read  sequentially  on  the  cards  in  a 
1015  format. 


22)  Plot  Variables 

A blank  card  must  be  included  here.  (This  option  when 
employed,  is  used  to  make  printer  plots  of  specified 
variables  when  a printer  plot  subroutine  is  available 
on  the  computer  system  being  used) . 

23)  Initial  Values  of  the  Unknown  Variables 


The  first  card  in  this  sequence  contains  two  values  in 
215  format.  These  numbers  (zero,  positive  or  negative) 
determine  the  manner  in  which  the  initial  values  of  the 
unknown  variables  and  their  derivatives  are  read.  (The 
first  number  pertains  to  the  variable  and  the  second 
number  pertains  to  the  derivative  of  the  variable). 

Zero 


If  a zero  is  read,  the  initial  value  of  all  variables 
is  read  in  a sequence  of  cards  in  5F10.9  format, 
hence,  the  number  of  cards  read  depends  upon  NOUT 
read  above. 

Positive  Number 


If  a positive  number,  say  N,  is  read  the  initial  value 
of  the  first  N variables  is  read  in  a sequence  of  cards 
in  5F10.9  format.  All  remaining  initial  values  are  then 
automatically  set  equal  to  zero. 


r 


i 


< 


Negative  Number 

If  a negative  number,  say  -N,  is  read,  the  initial 
values  of  N variables  are  read  randomly  on  N cards. 

Each  of  these  cards  contains  two  values.  The  first 
value  identifies  the  variable  number  and  the  second 
variable  is  the  initial  value  of  the  variable.  The 
format  for  these  two  values  is  15,  F15.9.  All  initial 
values  not  read  are  automatically  set  equal  to  zero. 

In  all  three  cases  above,  all  initial  values  of  the  variables  are 
read  before  the  initial  values  of  the  derivatives.  The  manner 
for  reading  the  initial  values  of  the  derivatives  is  the  same  as 
for  the  initial  values  of  the  variables.  (Mixed  modes  may  be 
employed) . 

The  usual  units  are  inches  and  inches/sec  for  displacement 
variables  and  degrees  and  degrees/sec  for  angular  variables. 

24)  Printing  Priority 

On  a single  card  in  15  format,  a number  is  read  between  -100 
and  +100  determining  the  amount  of  print-out  for  a given  run 
as  follows: 

If  0 all  print-out  (the  "standard  print-out")  for  every 
printing  time  (see  26  below)  is  given. 

If  > 10  printing  of  X and  DX  in  radians  is  omitted. 

If  > 20  printing  of  absolute  velocity  and  acceleration 
of  all  bodies  is  omitted. 

If  > 30  printing  of  position  of  all  bodies  is  omitted. 

If  > 40  printing  joint  moments  is  omitted. 

If  > 50  printing  of  DDX  in  degrees  is  omitted. 

If  > 95  All  printing  Is  omitted. 

If  < 0 (negative)  the  corresponding  print-out  is  given 
for  all  Integration  steps,  and  all  print-out  at  every 
print  time  (see  26  below)  is  given.  (For  example  -100 
gives  the  maximum  possible  print-out). 

25)  Printing  Time 

This  is  a value  read  on  a single  card  in  F10.9  format  determining 
the  time  increment  for  print-out.  This  may  be  the  same  or 
different  (usually  greater)  than  the  integration  time  increment 
PRMT(3) . 


26)  PRMT 
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SAMPLE  LISTING  OF  : 

INPUT  DATA 

1. 

* Error  Print  - - 

. 

2. 

*0.916 

0.916 

0.75 

0.178 

0. 

143-  - -i 

3. 

*0.178 

0.143 

0.433 

0.555 

0. 

358  j 

A. 

*0.555 

0.358-  - 

- 

- - 

_____  J 

5. 

*15.7 

0.0 

0.0 

6. 

*0.0 

9.593 

0.0 

I, 

7. 

*0.0 

0.0 

15.33 

1 

8. 

*15.7 

0.0 

0.0 

9. 

*0.0 

9.595 

0.0 

I, 

10. 

*0.0 

0.0 

15.33 

z 

11. 

*11.225 

0. 

0.0 

12. 

*0. 

7.74 

0. 

I, 

13. 

*0. 

0. 

11.163 

J 

14. 

*.0087 

15. 

* 

2.8156 

X4 

16. 

* 

2.8156 

17. 

*3.059 

18. 

*0. 

3.059 

0. 

ls 

19. 

*0.0 

0.0 

.0588 

5 * 

Inertia  Dyadics 

20. 

*.0087 

21. 

* 

2.8156 

22. 

* 

2.8156 

o 

23. 

*8.059 

24. 

*0. 

3.059 

0. 

I7 

25. 

*0.0 

0.0 

;0588 

26. 

*4.643 

0. 

0. 

27. 

*0. 

4.643 

0. 

X8 

28. 

*0. 

0. 

2.302 

29. 

*0.065 

0.0 

0.0 

30. 

*0. 

10.162 

0. 

Ip 

31. 

*0.0 

0.0 

10.162 

32. 

*12.2315 

J. 

0. 

• 

33. 

*0. 

12.388 

0. 

X10 

34. 

*0. 

0. 

0.1711 

35. 

*0.085 

0.0 

0.0 

36. 

*0. 

10.162 

0. 

X11 

37. 

*0.0 

0.0 

10.162 

38. 

*12.2316 

0. 

0. 

39. 

*0. 

12.388 

0. 

X12 

40. 

*0. 

0. 

0.1711 

Error  Set-Blank 

Masses  of  the 
12  Bodies 


*0. 

0. 

0.  - 

*0.0 

0.0 

4.05 

*0.0 

0.0 

3.95 

*4.46 

* 

7.025 

*4.46 

* 

7.025 

*0. 

0. 

6.35 

*9.8738 

0.0 

0.0 

*0. 

0. 

-9.0 

*9.8738 

0.0 

0.0 

*0. 

0. 

-9.0 

*0.0 

0.0 

4.05 

*0. 

0. 

8.1 

* 0. 

8.35 

5.8 

*11.7 

*0. 

-8.35 

5.8 

*11.7 
* 0. 

0. 

7.9 

* 0. 

3.07 

-.65 

*18.6 

0.0 

0.0 

* 0. 

-3.07 

-.65 

*18.6 

0.0 

0.0- 

R Vectors 


XI  Vectors 


12  Cockpit  Surfaces- 


65. 

*18.0 

12.0 

24.0 

1.0 

0.1 

66. 

*18.0 

0.0 

24.0 

1.0 

0.0 

67. 

*18.0 

-12.0 

24.0 

1.0 

-0.1 

68. 

*0.0 

18.0 

0.0 

0.0 

1.0 

69. 

*0.0 

18.0 

0.0 

0.0 

1.0 

70. 

*0.0 

-18.0 

0.0 

0.0 

-1.0 

71. 

*0.0 

-18.0 

0.0 

0.0 

-1.0 

72. 

*0.0 

0.0 

36.0 

0.0 

0.0 

73. 

*30.0 

0.0 

0.0 

1.0 

0.0 

74. 

*20.0 

0.0 

8.0 

0.70711 

0.0 

75. 

*20.0 

0.0 

8.0 

1.0 

0.0 

76. 

*20.0 

0.0 

4.0 

0.0 

0.0 

77. 

78. 

* No  Mag  Tape  - - 

♦Sample  Data  - - - 

- - -Magnetic  Tape  Card 

Heading  Card 

79. 

*2010.0 

18.6-  - 

. - - - - 

-Floor  Spring  And  Distance 

80. 

* 1 - - 

- - - 

. - - - - 

~ 

0.3 

0.3 

0.3 

0.0 

0.0 

0.0  Intrusion 
0.0  Surfaces 
1.0 
0.0 

-0.70711 

0.0 

1.0-  - -»  • 


82. 

*0.0 

0.0 

525.0 

0.0  4: 

83. 

*0.003 

-11592.0 

\lnitial 

84. 

*0.006 

-3864.0 

Velocity 

85. 

*0.01 

-15456.0 

86. 

*0.015 

-3864.0 

87. 

*0.02 

0.0 

88. 

*0.025 

-4636.8 

89. 

*0.07 

-5796.0 

90. 

*0.09 

-9660.0 

91. 

*0.1 

0.0  - - - 

- - . 

92. 

*600.0 

10.0-  - ■ 

- - - - -Neck  Spr 

93. 

*30.0 

70.0 

5.0 

-30.0 

94. 

*4.5 

4.5 

4.5 

-4.5 

95. 

*120.0 

85.0 

90.0 

0.0 

96. 

*0. 

90.0 

0. 

0. 

_Initial 

Displacenen 


Acceleration  Input  For 
The  Cockpit 


-90.0 

-4.5 

-60.0 

-45.0 


-5.0  - - - 

-4.5 
-45.0 
0.  Orientat'i 


r 


-120.0 


-60.0 

-45.0 


-45.0 

0. 


Angle 

Limits 


99. 

*45. 

80. 

60. 

-45. 

-60. 

-60. 

100. 

*15. 

45. 

45. 

-15. 

-80. 

-5. 

101. 

*0. 

30. 

0. 

0. 

-90. 

0. 

102. 

*15. 

45. 

5. 

-15. 

-80. 

-45. 

103. 

*0. 

30. 

0. 

0. 

-90. 

0.  - 

104. 

*10.0-  - - 

- - - 

- -Damping  Coefficient 

For 

The  Joints 

105. 

*2000000.  - 

- - - 

- -Damping  Coefficient 

For 

The  Joint 

Limit; 

106. 

*639.60 

639.60  - 

- -Spring  Constants  For  The  Seat 

107. 

*4.0 

4.0  - - 

- -Displacement  Limits 

For 

The  Seat 

108. 

* 2 Belts 

- - - 

- -Seat  Belts 

109. 

* 2 - - 

- - - 

- - 

- - 

- 

110. 

*4.0 

6.85 

3.0 

(Q) 

111. 

*-12. 

7.85 

-7. 

(RHO) 

112. 

*4296.0 

113. 

* 2 

Belt  Data 

1 

114. 

*4.0 

-6.85 

3.0 

115. 

-12. 

-7.85 

-7. 

116. 

*4296.0  - - 

- - - 

. _ - 

- - 

_ J 

1 Belt 

1 Seat 
1 Gravity 

1 Angle  Constraints 
No  Floor  Force 
No  Impulse 


Blank  Card 


Switches 


124. 

125. 

* 

* 

No  Neck 
1 Surf . - 

...J 

126. 

* 

28 

- - 

- - - 

- - - 

- - - 

- - Nout 

127. 

* 

7 8 

9 

13 

14 

15 

16  17  18 

19 

128. 

* 

20  21 

23 

25 

26 

27 

29  31  32 

33 

129. 

* 

34  35 

36 

38 

40 

41 

42  44 

- - 

JFXI  ARRAY 


0 No  Plots 

* -5  2 If  Read  Y.DY -t 

* 3 -1.84-^  | 

* 11  77.350  _ • „ . .. 

* 15  77.350  V Y mitlai  Y.lu.,  Of 

^ 22  ^ 360  i The  Uxilciiowii  Vctricibl.es 

* 26  -6.360 J 

*0.0  0.0  DY • 

* 11  Print-  --------------  Printing  Priority 

*0.005  Del  Print-  ----------  Printing  Time 

*0.0  2.00  .1953125-2.01  PRMT 

**End  Of  Data 


Blank  Card  For  Plots 


Initial  Values  Of 
The  Unknown  Variables 


OUTPUT  DATA 


The  output  data  is  labeled  on  the  computer  print-out  and  thus  is 
primarily  self-explanatory.  It  consists  of  two  parts:  The  first  part 

is  simply  an  "echo"  or  copy  of  the  input  data.  (In  the  input  motion 
section,  the  velocity  and  position  at  the  beginning  of  each  acceleration 
interval  is  computed  and  also  listed.)  The  second  part  contains  at  each 
output  time  (depending  upon  the  printing  priority  (See  24  and  25  in  the 
foregoing  section.): 

i)  All  variables  and  the  first  derivative  of  the 
variables  (angles  are  in  radians). 

ii)  Joint  positions  (in  both  inertia  space  and 
relative  to  the  vehicle). 

iii)  Mass  center  positions  (in  both  inertia  space  and 
relative  to  the  vehicle). 

iv)  Mass  center  velocities  and  accelerations. 

v)  All  variables  and  the  first  and  second 

derivatives  of  the  variables  (angles  are  in 
degrees) . 

vi)  Moments  and  forces  associated  with  variables 
which  are  specified. 

vii)  Seat  belt  forces. 

A sample  listing  of  output  data  at  a particular  time  interval  is  as 
follows : 
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Name:  UCIN-CRASH 

Developers:  C.E.  Passerello,  M.W.  Harlow,  and  R.L.  Huston,  Depart- 

ment of  Engineering  Analysis,  University  of  Cincinnati, 
Cincinnati,  Ohio. 

Date:  May  1,  1974 

Purpose  and  Use:  To  study  the  dynamics  of  crash  victims. 

Language:  FORTRAN  IV 

Operation:  Batch  or  time-sharing  (TYMSHARE) . 

Hardware:  Written  and  used  with  IBM  System  370. 

Documentation  and  Verification:  See  References  [15-18]. 

User  Manual:  See  Reference  [16] . 

Availability:  Available  from  authors  on  tape  at  cost  of  tape, 

handling,  and  mailing. 
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